Coulomb Effects in Few-Body Reactions 
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The method of screening and renormalization is used to include the Coulomb interaction between the charged 
particles in the momentum-space description of three- and four-body nuclear reactions. The necessity for the 
renormalization of the scattering amplitudes and the reliability of the method is demonstrated. The Coulomb 
effect on observables is discussed. 
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I. INTRODUCTION 

The inclusion of the long-range Coulomb interaction in the 
description of the three- and four-particle scattering is a chal- 
lenging task in theoretical few -body nuclear physics. The long 
range of the Coulomb potential prevents the direct application 
of the standard scattering theory. There is a number of sugges- 
tions how to overcome this difficulty; most of them are based 
on the configuration-space framework fl_l|-|4£| and are limited to 
energies below three-body breakup threshold, while the oth- 
ers O-Ba] have not matured yet into practical applications. Up 
to now only few approaches led to the results above three- 
body breakup threshold. Those are configuration-space calcu- 
lations for proton-deuteron (p-d) elastic scattering using the 
Kohn variational principle [9] and the screening and renormal- 
ization method in the framework of momentum-space integral 
equations ffTol — Tl~3tl : the latter approach will be discussed here 
in more details. Very recently p-d results above three-body 
breakup threshold were also obtained using modified Faddeev 
equation in configuration space together with the dumping of 
particular Coulomb contributions 01411 . 



II. METHOD OF SCREENING AND RENORMALIZATION 



In nature the Coulomb potential w c is always screened at 
large distances. The comparison of the data from typical nu- 
clear physics experiments and theoretical predictions with full 
Coulomb is meaningfull only if the full and screened Cou- 
lomb become physically indistinguishable. This was proved 
in Refs. Il5l ITdTI where the screening and renormalization 
method for the scattering of two charged particles was pro- 
posed. We base our treatment of the Coulomb interaction on 
that idea. 

The standard scattering theory is formally applicable to 
the screened Coulomb potential w R , i.e., the Lippmann- 
Schwinger equation yields the two-particle transition matrix 



t R - w R + w R g t R 



(1) 



where go is the free resolvent. In the configuration-space rep- 
resentation we choose the screened Coulomb potential as 



w R (r) = w c (r) e- ( ' /s) " 



(2) 



with R being the screening radius and n controlling the 
smoothness of the screening. It was proven in Ref. IU5I1 that 
in the limit of infinite screening radius R the on-shell screened 



Coulomb transition matrix (screened Coulomb scattering am- 
plitude) (p'felp) with p' = p, renormalized by the infinitely 
oscillating phase factor z R l (p) = e 2 ®"^', approaches the full 
Coulomb amplitude (p'|fclp) in general as a distribution, i.e., 

Mmz^OO J rfWMPMp) = j d 2 $(p'\t c \p)<p® (3) 

for any test function tp(p) with the properties given in Ref. 
O15II : in particular, it must vanish for forward scattering, i.e., 
<p(ji') = 0. Reference [15] uses partial wave expansion, which 
for the full Coulomb amplitude itself converges only as a dis- 
tribution, and it therefore is unable to make any conclusions 
on the possible pointwise convergence. However, as argued 
already in Ref. [15], the convergence of the renormalized 
screened Coulomb scattering amplitude to the full Coulomb 
amplitude in the sense of distributions is sufficient for the de- 
scription of physical observables. For a fixed final state obser- 
vation direction p' the cross section is determined not directly 
by the scattering amplitude, but by the outgoing wave packet 



<p f (p') = j d 3 p<p'|5|p)^(p) 



<rp<p'Mp>^(p)lp= 



p=p> 



(4) 
(5) 



that is related to the initial wave packet ^/(p) by the S -matrix 
or the scattering amplitude. Thus, in the step from the scatter- 
ing amplitude to the cross section one has to go through the 
conceptual exercise of averaging the scattering amplitude over 
the initial state physical wave packet being peaked around the 
experimental beam momentum p, . Thus, on the rhs of Eq. (01 
<p'|?clP/) itself is picked out in that average. In addition, the 
outgoing wave packet is never observed in the forward direc- 
tion, i.e., p' + p,; the necessary property y,(p') = can there- 
fore always be fulfilled by the sharpness of the initial wave 
packet. In the practical calculations ifTTl [l2tl the above av- 
eraging is carried out implicitly, replacing the renormalized 
screened Coulomb amplitude in the R — » 00 limit by the full 



one, i.e., 



lim z R \p){p'\t R \p) -> (p'lfdp)- 



«->co 



(6) 



Since z R l (p) is only a phase factor, the above relations indeed 
demonstrate that the physical observables become insensitive 
to screening provided it takes place at sufficiently large dis- 
tances R and, in the R — > 00 limit, coincide with the corre- 
sponding quantities referring to the full Coulomb. Further- 

more, renormalization by z R 2 (pi) in the R — > 00 limit relates 
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also the screened and full Coulomb wave functions lfl7ll . i.e., 
lim (1 + g t R )\ 9 )z^(p) = |^ +) (p)>. (7) 



A. Three-particle scattering 

The screening and renormalization method based on the 
above relations can be extended to more complicated systems, 
albeit with some limitations. The systems of two- and three- 
particles interacting via pairwise strong short-range and scree- 
ned Coulomb potentials, v a and w aR , a — 1, 2, 3, is considered 

in Refs. oi mm. 

Here we extend our treatment to a more 
general case where an irreducible three-body force 

3 

V (3 ) = £ u a (8) 

is present; it is decomposed into three terms u a . The full re- 
solvent 

G (R> = [E+ iO - Ho - Y/vy + u y + <V>] _1 . (9) 

r 

with //o being the three-particle kinetic energy operator and E 
the available energy may be decomposed into channel resol- 
vents 

Gf = (E + iO-H -v a -w aR y 1 , (10) 

and the multichannel three-particle transition operator Uj® 
according to 

G^=5 pa G^ + GflJ^Gf\ (11) 

all operators depend parametrically on the Coulomb screening 
radius R. The full multichannel transition matrix U^J for elas- 
tic and rearrangement scattering is calculated from the integral 
equation 11911 

7 

+ u a + Y J UyGod + T^GoW™ (12) 

y 

that is a generalisation of the Alt, Grassberger, and Sandhas 
(AGS) equation [20] in the presence of the three-body force. 
Go = (E + iO - i/o) 1 is me f ree resolvent, 6p a = 1 - 8a a , 
and the two-particle transition matrix is derived from the full 
channel interaction v a + w aR , i.e., 

T ( a R) = (v a + w aR ) + (v a + w aR )G T a R \ (13) 

The on-shell matrix elements (bptf\Ug a \b a q) are amplitudes 
(up to a factor) for elastic (J3 - a) and rearrangement (J3 + a) 
scattering. The channel states l^q) are the eigenstates of the 
corresponding channel Hamiltonian H a = Ho + v a with the 
energy eigenvalue E. \b a q) is a product of the bound state 
wave function \b a ) for the pair a and the plane wave with the 



relative particle-pair a momentum q; the dependence on the 
discrete quantum numbers is suppressed in our notation. 

In order to isolate the screened Coulomb contributions to 
the transition amplitude that diverge in the infinite R limit, 
we use a decomposition of the full resolvent into alternative 
channel resolvents 

G aR = (E + iO - Ho - v a - w aR - W%?)-\ (14) 

where is the screened Coulomb potential between the 
spectator particle a and the center of mass (cm.) of the re- 
maining pair. The same screening function is used for both 
Coulomb potentials w aR and W£™. The corresponding transi- 
tion matrix 

1 aR vy aR + vy aR 1 aR > KiJ > 

is a two-body operator and therefore its on-shell and half-shell 
behaviour in the limit R — » oo is given by Eqs. <[3j and ((T). It 
relates the channel resolvents as 

G aR = G {R) +G {R) T™GfK (16) 
Thus, the full resolvent can alternatively be decomposed into 
G {R) = 6 Pa G aR + G m U^G aR (17) 
= 8 pa G™ + (f8 pa T2&* 
+Gf [1 + 7*fGf ]8g>[l + Gf r^]Gf , (18) 

where the reduced transition operator U (R J(Z) may be calcu- 
lated through the integral equation 

Ufa = 5pa(G~ R + v a ) + u a + 5 j3a r W aR 

+ YjifyyVy + "y + fy/^G^E/^ (19) 
y 

which is driven by the strong two- and three-body poten- 
tials v a and u a and the potential of three-body nature 'Wan = 
Yjy(5ayW yR - 8 ay Wy^). This potential "W aR accounts for the 
difference between the direct Coulomb interaction w yR and the 
auxiliary one that takes place between the charged par- 
ticle and the cm. of the remaining pair. When calculated 
between on-shell screened Coulomb states, ui is of short- 

' pa 

range, even in the infinite R limit. Equation (JTTJ, together 
with Eq. ( fTTb . gives a relation between full and reduced three- 
particle transition operators, i.e., 

U% = + [1 + 7^Cf]0£[l + Gf T™1 (20) 

= 8 fSa T a ;£ + (Uf a ) -8 fla T a f). (21) 

Thus, the three-particle transition operator U R a has a long- 
range part SpaT^x whereas the remainder lAJ - 8p a T^ is 
a short-range operator that is externally distorted due to the 
screened Coulomb waves generated by [1 + G (R) T^]. On- 
shell, both parts do not have a proper limit as R — > oo but 
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the limit exists after renormalization by an appropriate phase 
factor, yielding the transition amplitude for full Coulomb 

^/jq'iagM) = <MMKc IM> 

+ Km [Z;j tf)(b0n'\(!J]g - 5 pa T^)\b a q)Z a l (q)]. (22) 

The first term on the right-hand side of Eq. d22l is known an- 
alytically lfl5ll ; it corresponds to the particle-pair a full Cou- 
lomb transition amplitude that results from the implicit renor- 
malization of according to Eq. ©. The R — > 00 limit for 
the remaining part (C/£ - 5p a T z a ™) of the multichannel tran- 
sition matrix is performed numerically; due to the short-range 
nature of this term, as demonstrated in Eq. ( 1201 . the conver- 
gence with the increasing screening radius R is fast and the 
limit is reached with sufficient accuracy at finite R; further- 
more, it can be calculated using the partial-wave expansion. 
We emphasize that Eq. (1221 is by no means an approximation 
since it is based on the obviously exact identity (fJTJ where the 
R — > 00 limit for each term exists and is calculated separately. 

The renormalization factor for R — > 00 is a diverging phase 
factor 

Z aR (q) = e- 2i °^"\ (23) 

where <f> a R(q), though independent of the particle-pair relative 
angular momentum / in the infinite R limit, may be realized 
by 

®aR(q) = - 1m(q), (24) 

with the diverging screened Coulomb phase shift T]iR(q) corre- 
sponding to standard boundary conditions and the proper Cou- 
lomb one <Ti(q) referring to the logarithmically distorted pro- 
per Coulomb boundary conditions. For the screened Coulomb 
potential of Eq. (f2| the infinite R limit of <& a R(q) is known an- 
alytically, 

<S> a n(q) = *W9)Pn (2qR) - C/n], (25) 

where C ~ 0.5772156649 is the Euler number and < K a {q) 
is the Coulomb parameter. The form of the renormalization 
phase <$ aR (q) to be used in the actual calculations with finite 
screening radii R is not unique, but the converged results show 
independence of the chosen form of <f> a R(q). 

For breakup observables we follow a very similar strat- 
egy where the starting point is the AGS -like equation for the 
breakup operator 

C = C ^X 7f >G„t/£> 

y 

+ u a + Y J "r G od + Tf>G )U™ (26) 

y 

and its relation to the full resolvent, i.e., 

G<*> = G q U™gM. (27) 
In the same spirit, we introduce auxiliary Coulomb resolvent 

G R = (E + iO -H -J] w r*y l - (28) 

y 



that keeps only the screened Coulomb interaction. The pro- 
per three-body Coulomb wave function and its relation to the 
three-body screened Coulomb wave function generated by Gr 
is, in general, unknown. This prevents the application of the 
screening and renormalization method to the reactions involv- 
ing three free charged particles (nucleons or nuclei) in the final 
state. 

However, in the system of two charged particles and a neu- 
tral one, only the channel y — p, corresponding to the corre- 
lated pair of charged particles, contributes to Gr which sim- 
plifies to 

G R = Gq + G T p rGo, (29) 
T p r = WpR + w p rGqT p r, (30) 

making channel p the most convenient choice for the de- 
scription of the final breakup state. Thus, for the purpose of 
breakup, a decomposition of the full resolvent, alternative to 
Eq. (|27) is 

G (R) ^G R U^G a R (31) 
= G [l + T pR G ]U^[l + G^T^Gf , (32) 

where the reduced breakup operator 0^ may be calculated 
through quadrature 

UqJ = GaR + + "a + £(<V + "y^y*^, (33) 

y 

from the corresponding reduced operator Uj^(Z) of elas- 
tic/rearrangement scattering. On-shell, the reduced operator 
Oi (Z) calculated between screened Coulomb distorted ini- 
tial and final states is of finite range, though the two contri- 
butions in Eq. (1331 have slightly different range properties as 
discussed in Ref. 11311 . The relation between the full and re- 
duced breakup operators is 

tC = (1 + T pR G )U™(l + Gf T™). (34) 

In the full breakup operator U^(Z) the external distortions 
show up in screened Coulomb waves generated by (1 + 
Ga T c " R °') in the initial state and by (1 + T p rGq) in the fi- 
nal state; both wave functions do not have proper limits as 
R — > oo. Therefore the full breakup transition amplitude in 
the case of the unscreened Coulomb potential is obtained via 
the renormalization of the on-shell breakup transition matrix 
Uq in the infinite R limit 

(p'q'I^IM) = ^a[z?(p')<p'q'l^l*aq>^i(9)], 

(35) 

where p' is the relative momentum between the charged par- 
ticles in the final state, q' the corresponding particle-pair rel- 
ative momentum, and 

ZR (p') = e -2i«(p')Hn(2pR)-C/n] 

the final-state renormalization factor. The limit in Eq. d35l 
has to be performed numerically, but, due to the short-range 
nature of the breakup operator, the convergence with the in- 
creasing screening radius R is fast and the limit is reached 
with sufficient accuracy at finite R. 



4 




FIG. 1. Screening function w R (r)/wc(r) as function of the distance r 
between charged particles for characteristic values of the parameter 
n in Eq. ([2j: n = 1 (green curve) corresponds to Yukawa screening, 
n = 4 (red curve) is our standard choice, and n — » oo (blue curve) 
corresponds to a sharp cutoff. 
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FIG. 2. Convergence of the 1 Sq pp phase shift rj with screening 
radius R for 3 MeV proton lab energy. Results obtained with n = 1 
(green curve), n = 4 (red curve), and n — > oo (blue curve) screening 
functions are compared with exact value given by black line. 



B. Practical realization 

To calculate the short-range part of the elastic, rearrange- 
ment, and breakup scattering amplitudes d22l and d35l ) we 
solve standard scattering equations fl% , JTBt . and d26*i > at fi- 
nite Coulomb screening radius R using the momentum-space 
partial-wave representation as described in detail in Refs. 112 ll - 
12311 . We have to make sure that R is large enough to achieve 
(after renormalization) the /^-independence of the results up 



to a desired accuracy. However, those R values are larger 
than the range of the nuclear interaction resulting in a slower 
convergence of the partial-wave expansion. As we found in 
Ref. B12I1 . the practical success of the screening and renormal- 
ization method depends strongly on the choice of the scree- 
ning function, i.e., on the choice of the exponent n in Eq. (12). 
One of the essential differences compared to previous works 
ifioi [Til l24l l25ll is that we use a sharper screening than the 
Yukawa screening (n = 1). We want to ensure that the scree- 
ned Coulomb potential w R approximates well the true Cou- 
lomb one w c for distances r < R and simultaneously van- 
ishes rapidly for r > R, providing a comparatively fast con- 
vergence of the partial-wave expansion. However, the sharp 
cutoff (« — » 00) yields an unpleasant oscillatory behavior in 
the momentum-space representation, leading to convergence 
problems. We find values 3 < n < 8 to provide a sufficiently 
smooth, but at the same time a sufficiently rapid screening 
around r — R. The screening functions for different n val- 
ues are compared in Fig. Q] showing that the present choice 
n — 4 includes much more of the exact Coulomb potential 
at short distances than the Yukawa screening used previously. 
For example, Yukawa screening requires a screening radius of 
R = 1280 fm in order to approximate true Coulomb at rel- 
ative distance r — 5 fm as well as the present choice does 
with R — 20 fm. As shown in Fig. [2] for the proton-proton 
(pp) 1 So phase shift, the convergence with the screening ra- 
dius for n — 4 is much faster than for Yukawa screening or 
sharp cutoff. Furthermore, due to the much shorter range of 
the n = 4 screening compared to the Yukawa screening, the 
quasisingularities of the screened Coulomb potential are far 
less pronounced in the n — 4 case. With our optimal choice 
3 < n < 8 the convergence of the partial-wave expansion, 
though being slower than for the nuclear interaction alone, 
can still be achieved; alternatively, a perturbative approach for 
higher two-particle partial waves B26I1 that is computationally 
less demanding but, nevertheless, highly reliable as proved in 
Ref. 12711 . can be used as well. 

C. Renormalization in proton-deuteron scattering 

For the nucleon-deuteron scattering it is convenient to con- 
sider nucleons as identical particles and use the isospin for- 
malism. Since the isospin conservation is violated by the Cou- 
lomb force, both total 3A^ isospin T — I and | states have to 
be included. The symmetrized transition operator for elas- 
tic scattering is the solution of the symmetrized AGS integral 
equation lfl9ll 

U (R) = PGo 1 +(1+P)u + PT (R) G U iR) 

+(1 + P)uG (\ + T (R) G Q )U (R) . (37) 

We omit the spectator index that is not needed anymore. The 
basis states are antisymmetric in the pair only, the full anti- 
symmetry is ensured by P = PnPn + PnPih where P a p is 
the permutation operator of particles a and f3. The breakup 
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FIG. 3. Real and imaginary parts of nonrenormalized (left side) and renormalized (right side) p-d elastic scattering amplitudes (in arbitrary 
units) at E p = 10 MeV shown as functions of the cm. scattering angle. The initial and final particle spin projection quantum numbers are 
nip = 5, m' = —\, and m d = m' d = 1. Results for the CD Bonn + A two-baryon potential obtained with the screening radius R = 5 fm (black 
curves), 10 fm (blue curves), 15 fm (green curves), and 20 fm (red curves) are compared. 



operator is then obtained from the quadrature 

U (R> = (1 + p)[G ' + u + T (R) G U (R) + uG (l + T (R) G )U (R) ]. 

(38) 

The amplitudes for p-d elastic scattering and breakup refer- 
ring to unscreened Coulomb are obtained after the renormal- 
ization of the corresponding on-shell matrix elements of the 
symmetrized transition operators, i.e., 

<q'|f7 (C) lq) = (q'lrriQ) 

+ lim [Z^(q){q'\(U (R) ~ T c R m )\q)], (39) 

R— >oo 

<p'q'|t/< C) |q> = l^[zlhp'Wq'\U^\)Z~Hq)l (40) 

However, a recent work j28ll on an alternative Coulomb 
treatment in pd scattering proposed a different renormaliza- 
tion prescription: p-d elastic scattering amplitude calculated 
with screened Coulomb does not need renormalization at all, 
i.e., the limit limR^oo(q'|(7 (R) |q) should exist; the renormal- 



ization for the breakup amplitude is needed but it is differ- 
ent from ours given in Eq. ( |40l . However, numerical results 
of Ref. II28I1 involve approximations that are not well under 
control. First, the screened Coulomb proton-proton transition 
matrix is approximated by the screened Coulomb potential 
thereby replacing a complex quantity by a real one. Second, 
the contributions 
(q'\PT R PG T m GoU^\q) and 

<p'q'|(l + P)T R PGoT w G U {R) \q) to the elastic scattering and 
breakup amplitudes, respectively, are neglected; this may have 
even more serious consequences since the above terms among 
others include screened Coulomb contributions that are of the 
leading (first) order. Indeed, under the approximations of 
Ref. 12811 the ^-dependence of the resulting screened Coulomb 
amplitude for p-d elastic scattering is weak and the R — > 00 
limit of renormalized amplitude in Eq. ( f39l > does not exists. 
However, the situation changes completely when the scat- 
tering amplitudes are calculated exactly as in our work. In 
Figs. [3] and |4] we study the dependence on the screening ra- 
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FIG. 4. Real and imaginary parts of nonrenormalized (left side) and renormalized (right side) p-d elastic scattering amplitudes (in arbitrary 
units) at E p = 10 MeV shown as functions of the cm. scattering angle. The initial and final particle spin projection quantum numbers are 
m p = m' p = j, and m d = m' d = 1. Curves as in Fig. [3] 



dius R for selected components of the nonrenormalized and 
renormalized p-d elastic scattering amplitudes (q'|£/ (S) |q) and 
<q'|t/ (C) lq>- The CD Bonn + A two-baryon potential JH 
is taken as the hadronic interaction and the proton lab en- 
ergy E p = 10 MeV. In Fig. [3] we show spin-nondiagonal 
amplitudes that have only the Coulomb-distorted short-range 
part. The nonrenormalized amplitude (q'|t/ (S) |q) shows a 
clear /^-dependence; in fact, its absolute value becomes R- 
independent but the phase changes like InR. In contrast, the 
renormalized amplitude (q'|£/ (C) |q), within the accuracy of the 
plot, becomes independent of R for R > 10 fm. The spin- 
diagonal amplitude having also a long-range part shows, as 
expected, very strong /^-dependence at small scattering an- 
gles before renormalization but becomes /^-independent after 
renormalization as well as Fig. ^demonstrates. In Fig. |5]we 
show the corresponding results for p-d breakup at E p = 13 
MeV in the space star configuration. For breakup reaction 
the renormalized amplitude (p'n'\U^\q) converges with the 
screening radius somehow slower than for elastic scatter- 



ing, it becomes independent of R for R > 20 fm whereas 
the nonrenormalized amplitude (p'q'lt/tf'lQ) shows a clear 
/^-dependence. Furthermore, although the breakup ampli- 
tude (without the three-nucleon force) can be decomposed 
into three terms Ti ^GqUjJ, separately none of them has an 
R — > oo limit with or without renormalization in contrast to the 
conjecture of Ref. 12811 : infinite R limit exists only for renor- 
malized full breakup amplitude. Thus, fully converged nu- 
merical results without uncontrolled approximations clearly 
support the standard screening and renormalization theory as 
given in Eqs. (l39l > and (1401 and not the one of Ref. [28]. 

The convergence of the observables with the screening ra- 
dius R used to calculate the Coulomb-distorted short-range 
part of the amplitudes is the internal criterion for the relia- 
bility of our method. Numerous examples can be found in 
Refs. lfl2l [l3l [30l l3lll . In most cases the convergence is 
impressively fast; the screening radius R = 10 to 30 fm is 
sufficient. The exceptions requiring larger screening radii are 
the observables at very low energies and the breakup differ- 
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FIG. 5. Real and imaginary parts of nonrenormalized (left side) and renormalized by exp{i[K(p) + 'Kiq)] \nR] (right side) p-d breakup ampli- 
tudes (in arbitrary units) at E„ = 13 MeV in the space star configuration shown as functions of the arclength S along the kinematical curve. 
Curves as in Fig. [3] 



ential cross section in kinematical situations characterized by 
very low relative energy E K \ between the two charged parti- 
cles, e.g., p-d breakup or photodisintegration of 3 He close to 
the pp final-state interaction (pp-FSI) regime [13]. The slow 
convergence under those conditions is not surprising, since the 
renormalization factor itself as well as the Coulomb parameter 
become ill-defined, indicating that the screening and renor- 
malization procedure cannot be applied at E K \ = 0. There- 
fore an extrapolation has to be used to calculate the observ- 
ables at E re \ = 0, which works pretty well when the observ- 
ables vary smoothly with E K \ which is the case in three-body 
breakup reactions. Furthermore, the results for p-d elastic 
scattering obtained by the present technique were compared 
l32tl with those of Ref. H obtained from the variational solu- 
tion of the three-nucleon Schrodinger equation in configura- 
tion space with the inclusion of an unscreened Coulomb po- 
tential between the protons and imposing the proper Coulomb 
boundary conditions explicitly. Good agreement over a wide 
energy range was found indicating that both techniques for 



including the Coulomb interaction are reliable. At very low 
energies the coordinate-space treatments remain favored since 
there the method of screening and renormalization converges 
slowly and therefore becomes technically too demanding, but 
at higher energies and for three-body breakup reactions it is 
more efficient. 



HI. RESULTS 

In this section we present selected results for reactions in 
various three- and four-body nuclear systems. 

A. Proton-deuteron scattering 

Realistic models for the nuclear interaction are used, e.g., 
CD Bonn +A two-baryon potential H29I1 yielding an effec- 
tive three-nucleon force or the AVI 8 two-nucleon potential 



8 




60 120 180 60 120 180 60 120 180 

©cm. (deg) ©c. m . (cleg) ® c m . (deg) 



FIG. 6. Differential cross section and analyzing powers for p-d elastic scattering at 5 MeV proton lab energy as functions of the c.m. scattering 
angle. Results including the Coulomb interaction (red curves) are compared to results without Coulomb (blue curves). Hadronic potential is 
CD Bonn + A. The experimental data are from Ref. l33ll . 




FIG. 7. Differential cross section for p-d breakup at 13 MeV proton lab energy in space star (left), quasifree scattering (middle), and collinear 
(right) configurations as function of the arclength S along the kinematical curve. Results including the Coulomb interaction (red curves) are 
compared to results without Coulomb (blue curves). Hadronic interaction model is AVI 8 + UIX. The experimental data are from Ref. l34ll 
(full circles). For the space star configuration also the nd data from Refs. l35Ll36ll (open and full squares) are shown. 



1 37] with the irreducible Urbana IX (UIX) three-nucleon force deuteron vector analyzing powers. In p-d elastic scattering at 

1 3811 . The Coulomb effect is important in low-energy p-d elas- higher energies the Coulomb effect gets confined to forward 

tic scattering as Fig. ^demonstrates; the description of the ex- angles. 

perimental data is quite satisfactory except for the proton and In Fig. [7] we show results for p-d breakup at 13 MeV pro- 
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ton lab energy. Although the inclusion of the Coulomb force 
slightly improves the agreement with data in the space star 
configuration, the Coulomb effect is far too small to reproduce 
the difference between the experimental p-d and n-d data and 
to resolve the so-called space star anomaly. Slightly larger 
and beneficial Coulomb effects are seen in quasifree scatter- 
ing (QFS) and collinear configurations. 

The Coulomb effect in may become large even at higher 
energies when the relative p p energy in the final state is very 
low as demonstrated in Fig. [8] for p-d breakup at 130 MeV 
deuteron lab energy. In there, the Coulomb repulsion is re- 
sponsible for decreasing the cross section; the pp-FSl peak 
obtained in the absence of Coulomb may even be converted 
into a minimum. However, in some other configurations Cou- 
lomb may lead to a moderate increase of the differential cross 
section as shown in Fig. [9] 

B. Photodisintegration of 3 He 

The results are obtained with the CD Bonn +A two-baryon 
potential and include effective two- and three-nucleon elec- 
tromagnetic currents mediated by the A isobar l4ll l42ll . As 
in p-d breakup, the Coulomb effect may be very strong in 
the three-body photodisintegration of 3 He close to the pp-FSl 
kinematics. The effect remains also after partial integration 
when calculating the semi-inclusive observables as shown in 
Figs. [TOl and [TTI Not only the differential cross section gets 
affected but also the spin observables like the linear photon 
asymmetry for 3 He(y, p)pn reaction. 

C. Three-body nuclear reactions 

The screening and renormalization method in the frame- 
work of momentum-space AGS equations can be applied to 
the description of three-body-like nuclear reactions involving 
heavier nuclei whose interaction with nucleons is described by 
the optical potentials. Examples are deuteron scattering on a 
stable nucleus A or proton scattering on a weakly bound two- 
body system (An) consisting of the core A and the neutron n; 
all elastic, transfer, charge-exchange, and breakup reactions, 
allowed by the chosen Hamiltonian, can be calculated on the 
same footing. The calculations have been performed for a 
number of reactions with the core A ranging from the a parti- 
cle to 58 Ni [30, HQ. 

In Fig. [12] we show the differential cross section for a-d 
breakup reaction at 15 MeV a lab energy that is calculated us- 
ing three-body (a, p, n) model. The most important Coulomb 
effect is the shift of the ap P-wave resonance position that 
leads to the corresponding changes in the structure of the ob- 
servables. The predictions without Coulomb fail completely 
in accounting for the experimental data, while inclusion of 
the Coulomb moves the peaks of the differential cross sec- 
tion to the right positions, although the height of those peaks 
is not always reproduced, probably due to deficiencies of the 
nucleon-a potentials. 



The developed technique allowed to test the accuracy of 
traditional approximate nuclear reaction approaches like Con- 
tinuum Discretized Coupled Channels (CDCC) method ll47tl 
as has been done in Ref. Q48I] for d + 12 C, d + 58 Ni, and 
p + n Be reactions. It was found that CDCC is indeed a re- 
liable method to calculate d + A elastic and breakup cross 
sections but may lack accuracy for transfer reactions such as 
p + n Be — > d + I0 Be and for breakup of one-neutron halo nu- 
clei p + n Be — » p + n + 10 Be. Furthermore, novel dynamical 
input like energy-dependent M49I1 or nonlocal optical poten- 
tials ll50l BTll could be included for the first time due to the 
use of the momentum-space framework. Especially important 
nonlocality effects were found for (d, p) and (p, d) transfer re- 
actions involving stable [52[] as well as exotic nuclei lf53ll . 

D. Four-nucleon scattering 

Exact description of the four-nucleon scattering is given by 
the Faddeev-Yakubovsky equations [54] for the wave-function 
components or by the equivalent AGS equations JUt]. We use 
isospin formalism and solve the symmetrized AGS equations 
ll56ll for the transition operators. Coulomb interaction is in- 
cluded using the same idea of screening and renormalization. 
Long- and Coulomb-distorted short-range parts in the scatter- 
ing amplitudes are separated 1571 l5al . The former is of two- 
body nature and its R — > 00 limit is known analytically. The 
Coulomb-distorted short-range part is calculated by solving 
symmetrized AGS equations numerically where the screened 
Coulomb potential is added to the nuclear proton-proton po- 
tential; the R — > 00 limit is reached with sufficient accuracy at 
finite screening radii R as demonstrated in Ref. f57ll . 

The four-nucleon scattering calculations have been per- 
formed so far only below three-body breakup threshold where 
the Coulomb effect is extremely important. This is demon- 
strated in Fig. [14] for proton- 3 He elastic scattering at 4 MeV 
proton lab energy. As in p-d scattering, the proton analyzing 
power is underpredicted by the theory i57[ 15911 . 

In Figs. Q3] and [16] we show examples for transfer and 
charge-exchange reactions in the four-nucleon system. The 
two-nucleon interactions we use are AVI 8 Q37[| , the one de- 
rived from chiral perturbation theory at next-to-next-to-next- 
to-leading order (N3LO) JH, CD Bonn (H, and inside non- 
local outside Yukawa (INOY04) potential by Doleschall JH]. 
Although here we do not include a three-nucleon force, its 
presence is simulated by using the potential INOY04 that fits 
both 3 He and 3 H experimental binding energies (7.72 MeV 
and 8.48 MeV, respectively). The results for the two-baryon 
potential CD Bonn + A M29I1 allowing for a virtual excitation 
of a nucleon to a A-isobar and thereby yielding consistent ef- 
fective three- and four-nucleon forces are qualitatively simi- 
lar and can be found in Ref. 06511 . Most of the experimental 
data are quite well described ll57l - l59[ l66ll at least by some of 
the used two-nucleon force models, but there exist also sev- 
eral discrepancies, e.g., for the neutron- 3 H total cross section 
ll56il67ll and for the proton analyzing power in the proton- 3 He 
elastic scattering J53, HI, H and in the p + 3 H -h> n + 3 He 
charge-exchange reaction [58]. 
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FIG. 8. Differential cross section for pd breakup at 130 MeV deuteron lab energy in selected kinematical configurations with small relative 
azimuthal angle. Curves as in Fig.|7]and the experimental data from Refs. 1 39, 40], 
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FIG. 9. Differential cross section for p-d breakup at 130 MeV deuteron lab energy in selected kinematical configurations with large relative 
azimuthal angle. Curves as in Fig.[7]and the experimental data from Ref. 1 39J]. 



IV. SUMMARY 

We used the method of screening and renormalization to in- 
clude the Coulomb interaction between the charged particles 



in few-body nuclear reactions. We demonstrated analytically 
and numerically that the limit of infinite screening radius ex- 
ists only for renormalized amplitudes. The short range part of 
the scattering amplitudes were obtained from the exact few- 
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FIG. 11. The semi-inclusive linear photon asymmetry for the 3 He(y, n)pp (left) and 3 He(y, p)pn (right) reactions at 15 MeV photon lab energy 
as function of the detected nucleon energy at 90° nucleon scattering angle. Curves as in Fig. [6] 
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FIG. 12. Differential cross section for a-d breakup at 15 MeV a lab energy in selected kinematical configurations as function of the arclength 
S along the kinematical curve. Results including the Coulomb interaction (red curves) are compared to results without Coulomb (blue curves). 
The experimental data are from Ref. 14311 . 



body scattering equations that were solved in the momentum- three- and four-nucleon scattering and for three-body-like nu- 
space framework. We obtained fully converged results for clear reactions. 



[1] A. Kievsky, S. Rosati, W. Tornow, M. Viviani, Nucl. Phys. 

A607, 402 (1996) 
[2] C.R. Chen, J.L. Friar, G.L. Payne, Few-Body Syst. 31, 13 

(2001) 

[3] S. Ishikawa, Few-Body Syst. 32, 229 (2003) 
[4] P. Doleschall, Z. Papp, Phys. Rev. C72, 044003 (2005) 
[5] E.O. Alt, S.B. Levin, S.L. Yakovlev, Phys. Rev. C69, 034002 
(2004) 

[6] A.S. Kadyrov, I. Bray, A.M. Mukhamedzhanov, A.T. Stel- 

bovics, Phys. Rev. A 72, 032712 (2005) 
[7] A.S. Kadyrov, I. Bray, A.M. Mukhamedzhanov, A.T. Stel- 

bovics, Ann. Phys. 324, 1516 (2005) 
[8] S. Oryu, Phys. Rev. C73, 054001 (2006) 
[9] A. Kievsky, M. Viviani, S. Rosati, Phys. Rev. C64, 024002 

(2001) 



[10] E.O. Alt, M. Rauh, Few-Body Syst. 17, 121 (1994) 

[11] E.O. Alt, A.M. Mukhamedzhanov, M.M. Nishonov, A.I. Sat- 

tarov, Phys. Rev. C 65, 064613 (2002) 
[12] A. Deltuva, A.C. Fonseca, PU. Sauer, Phys. Rev. C71, 054005 

(2005) 

[13] A. Deltuva, A.C. Fonseca, PU. Sauer, Phys. Rev. C72, 054004 
(2005) 

[14] S. Ishikawa, Phys. Rev. C80, 054002 (2009) 

[15] J.R. Taylor, Nuovo Cimento B 23, 313 (1974) 

[16] M.D. Semon, J.R. Taylor, Nuovo Cimento A 26, 48 (1975) 

[17] V.G. Gorshkov, Sov. Phys.-JETP 13, 1037 (1961) 

[18] A. Deltuva, A.C. Fonseca, PU. Sauer, 

Annu. Rev. Nucl. Part. Sci. 58, 27 (2008) 
[19] A. Deltuva, Phys. Rev. C80, 064002 (2009) 



13 




60 120 180 60 120 180 



©cm. (deg) © c . m . (cleg) 

FIG. 13. Differential cross section for d + 16 — » p + n O transfer to the 17 ground state 5/2 + (left side) and excited state 1 /2 + (right side) at 
Ei = 25.4 and 36.0 MeV calculated with local and nonlocal optical potentials. The experimental data are from Ref. l44Tl . 




FIG. 14. Differential cross section and proton analyzing power for p- 3 He elastic scattering at 4 MeV proton lab energy as functions of the cm. 
scattering angle. Results including the Coulomb interaction (red curves) are compared to results without Coulomb (blue curves). Hadronic 
potential is CD Bonn. The experimental data are from Refs. I60l.l6lll . 



[20] E.O. Alt, P. Grassberger, W. Sandhas, Nucl. Phys. B2, 167 
(1967) 

[21] K. Chmielewski, A. Deltuva, A.C. Fonseca, S. Nemoto, PU. 

Sauer, Phys. Rev. C 67, 014002 (2003) 
[22] A. Deltuva, K. Chmielewski, PU. Sauer, Phys. Rev. C67, 

034001 (2003) 



[23] A. Deltuva, Ph.D. thesis, University of Hannover (2003), 

|http://edokO 1 .tib.um-hannover~d e/edoks/eO 1 dh0 3/374454701.pdf| 
[24] GH. Berthold, A. Stadler, H. Zankel, Phys. Rev. C41, 1365 
(1990) 

[25] E.O. Alt, A.M. Mukhamedzhanov, A.I. Sattarov, 
Phys. Rev. Lett. 81, 4820 (1998) 



14 




50 100 150 50 100 150 

c.m. ( de 9) c.m. ( de 9) 



FIG. 15. Differential cross section and proton analyzing power of p + 3 H — » n + 3 He reaction at 2.48 and 6 MeV proton lab energy calculated 
with various realistic two-nucleon potentials. The cross section data are from Refs. 1 68] (circles) and 1 69] (squares) at 2.48 MeV, and from 
Ref. 1 73] at 6 MeV. A y data are from Ref. Hfl] at 2.48 MeV and from Ref. 1 72] at 6 MeV. 



[26] A. Deltuva, K. Chmielewski, P.U. Sauer, Phys. Rev. C67, 
054004 (2003) 

[27] A. Deltuva, A.C. Fonseca, P.U. Sauer, Phys. Rev. C73, 057001 
(2006) 

[28] H. Witala, R. Skibinski, J. Golak, W. Glockle, Eur. Phys. J. 

A41, 369 (2009); A41, 385 (2009) 
[29] A. Deltuva, R. Machleidt, P.U. Sauer, Phys. Rev. C68, 024005 

(2003) 

[30] A. Deltuva, Phys. Rev. C74, 064001 (2006) 
[31] A. Deltuva, A.C. Fonseca, P.U. Sauer, Nucl. Phys. A790, 52c 
(2007) 

[32] A. Deltuva, A.C. Fonseca, A. Kievsky, S. Rosati, P.U. Sauer, 

M. Viviani, Phys. Rev. C71, 064003 (2005) 
[33] K. Sagara, H. Oguri, S. Shimizu, K. Maeda, H. Nakamura, 

T. Nakashima, S. Morinobu, Phys. Rev. C50, 576 (1994) 
[34] G. Rauprich, S. Lemaitre, P. Niessen, K.R. Nyga, R. Reck- 

enfelderbaumer, L. Sydow, H. Paetz gen. Schieck, H. Witala, 

W. Glockle, Nucl. Phys. A535, 313 (1991) 
[35] J. Strate et al, Nucl. Phys. A501, 51 (1989) 
[36] H.R. Setze et al., Phys. Rev. C71, 034006 (2005) 
[37] R.B. Wiringa, V.G.J. Stoks, R. Schiavilla, Phys. Rev. C51, 38 

(1995) 

[38] B.S. Pudliner, VR. Pandharipande, J. Carlson, S.C. Pieper, R.B. 

Wiringa, Phys. Rev. C56, 1720 (1997) 
[39] S. Kistryn et al, Phys. Rev. C72, 044006 (2005) 



[40] S. Kistryn et al, Phys. Lett. B 641, 23 (2006) 

[41] A. Deltuva, L.P. Yuan, J. Adam Jr., A.C. Fonseca, P.U. Sauer, 

Phys. Rev. C 69, 034004 (2004) 
[42] A. Deltuva, L.P. Yuan, J. Adam Jr., P.U. Sauer, Phys. Rev. C70, 

034004 (2004) 

[43] I. Koersner, L. Glantz, A. Johansson, B. Sundqvist, H. Naka- 
mura, H. Noya, Nucl. Phys. A286, 431 (1977) 

[44] M.D. Cooper, W.F. Hornyak, PG. Roos, Nucl. Phys. A218, 249 
(1974) 

[45] A. Deltuva, Nucl. Phys. A821, 72 (2009) 
[46] E. Cravo, R. Crespo, A. Deltuva, A.C. Fonseca, Phys. Rev. C79, 
064610 (2009) 

[47] N. Austern, Y. Iseri, M. Kamimura, M. Kawai, G. Rawitscher, 

M. Yahiro, Phys. Rep. 154, 125 (1987) 
[48] A. Deltuva, A.M. Moro, E. Cravo, F.M. Nunes, A.C. Fonseca, 

Phys. Rev. C 76, 064602 (2007) 
[49] A. Deltuva, A.C. Fonseca, Phys. Rev. C79, 014606 (2009) 
[50] M.M. Giannini, G. Ricco, Ann. Phys. (NY) 102, 458 (1976) 
[51] M.M. Giannini, G. Ricco, A. Zucchiatti, Ann. Phys. (NY) 124, 

208 (1980) 

[52] A. Deltuva, Phys. Rev. C79, 021602 (2009) 
[53] A. Deltuva, Phys. Rev. C79, 054603 (2009) 
[54] O.A. Yakubovsky, Yad. Fiz. 5, 1312 (1967) [Sov. J. Nucl. Phys. 
5, 937 (1967)] 




FIG. 16. Differential cross section and deuteron analyzing powers of d + d — > p + 3 H and d + d —* n + 3 He reactions at 3 MeV deuteron lab 
energy. The cross section data are from Refs. l73ll (squares) and 1 74] (circles) and the analyzing power data are from Refs. (74l . l75ll . 



16 



[55] P. Grassberger, W. Sandhas, Nucl. Phys. B2, 181 (1967); E. O. 
Alt, P. Grassberger, and W. Sandhas, JINR report No. E4-6688 
(1972) 

[56] A. Deltuva, A.C. Fonseca, Phys. Rev. C75, 014005 (2007) 
[57] A. Deltuva, A.C. Fonseca, Phys. Rev. Lett. 98, 162502 (2007) 
[58] A. Deltuva, A.C. Fonseca, Phys. Rev. C76, 021001 (2007) 
[59] M. Viviani, A. Kievsky, S. Rosati, E.A. George, L.D. Knutson, 

Phys. Rev. Lett. 86, 3739 (2001) 
[60] B.M. Fisher, C.R. Brune, H.J. Karwowski, D.S. Leonard, E.J. 

Ludwig, T.C. Black, M. Viviani, A. Kievsky, S. Rosati, Phys. 

Rev. C74, 034001 (2006) 
[61] M.T. Alley, L.D. Knutson, Phys. Rev. C48, 1890 (1993) 
[62] D.R. Entem, R. Machleidt, Phys. Rev. C68, 041001(R) (2003) 
[63] R. Machleidt, Phys. Rev. C63, 024001 (2001) 
[64] P. Doleschall, Phys. Rev. C69, 054001 (2004) 
[65] A. Deltuva, A.C. Fonseca, P.U. Sauer, Phys. Lett. B660, 471 

(2008) 

[66] R. Lazauskas, Phys. Rev. C79, 054007 (2009) 



[67] R. Lazauskas, J. Carbonell, Phys. Rev. C70, 044002 (2004) 
[68] M. Drosg, Nucl. Sci. Eng. 67, 190 (1978) 
[69] G.A. Jarvis, Los Alamos Scientific Lab Reports No. 2014, 35 
(1956) 

[70] W.E. Wilson, R.L. Walter, D.B. Fossan, Nucl. Phys. 27, 421 
(1961) 

[71] M.A. Doyle, H.W. Clark, L.J. Dries, J.L. Regner, T.R. 

Donoghue, G.M. Hale, Nucl. Phys. A371, 225 (1981) 
[72] J.J. Jarmer, R.C. Haight, J.E. Simmons, J.C. Martin, T.R. 

Donoghue, Phys. Rev. C9, 1292 (1974) 
[73] J.M. Blair, G. Freier, E. Lampi, W. Sleator, J.H. Williams, Phys. 

Rev. 74, 1599 (1948) 
[74] W. Gruebler, V. Konig, PA. Schmelzbach, R. Risler, R.E. 

White, P. Marmier, Nucl. Phys. A193, 129 (1972) 
[75] L.J. Dries, H.W. Clark, R. Detomo, T.R. Donoghue, Phys. Lett. 

80B, 176 (1979) 



